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Method and Device for Gas Lifted Wells 
Field of invention 

An automatic model-based controller is used for stabili- 
5 zation of gas lifted oil wells. The controller stabilizes 
the pressures, temperatures and flow rates in the well 
through active feedback control and continuous manipula- 
tion of the opening of the production choke and/or the 
opening of the gas injection choke as a dynamic function 
10 of available process measurements. 

The system comprises a production choke for control of 
the product from the oil well and/or a gas injection 
choke for control of lift-gas to the annulus . 

15 

Background 

Unstable production conditions often occur for oil wells 
where gas lift is used to increase the oil production. 
This is a serious problem due to the fact that the insta- 
20 bilities often occur in the most optimal area of opera- 
tion for the gas lifted wells, thus causing the well to 
produce less than the system design capacity. 

Unstable production of gas lifted wells is not in agree - 
25 ment with smooth operation and it implies safety aspects 
and shutdown risks. The total oil and gas production must 
usually be less than the systems design capacity (e.g. of 
the separators) to allow for the peak production. Unsta- 
ble operation decreases sharply the lift gas efficiency, 
3 0 and leads to difficulties with gas lift allocation compu- 
tation . 



WO 00/00715 PCT/NO99/00217 ' 

2 

Unstable production of gas lifted wells may be caused by 
a large variety of factors, such as incorrect gas lift 
string design, improper valve setting, wrongly sized in- 
jection valve port, variation in supply pressure, or 
5 valve leaking or plugging. It is often difficult to find 
the origin of the instabilities. As a result, a pragmatic 
approach has often been used to solve the problem of un- 
stable production in short term. For example, if unstable 
production occurs, the operator often increases the 

10 amount of lift gas or increases the back-pressure by ad- 
justing the wellhead production choke to a smaller open- 
ing (choking) . Although these methods can be effective in 
reducing the instabilities, the production is still inef- 
ficient as either too much lift gas is used (high cost 

15 and limited availability of lift gas) or the well is pro- 
ducing against a high back pressure (at low rate) . In 
most cases, too much gas is injected into the gas lifted 
wells or the production rate is not maximized. 

2 0 Summary of the Invention 

The present invention relates to a method and a stabiliz- 
ing well controller for stabilization of gas lifted oil 
wells without using lift-gas in an inefficient way or by 
introducing a high static back pressure. The characteriz- 
25 ing features of the method are given in claim 1. The 

characterizing features of the stabilizing well control- 
ler are given in claim 14 . 

This concept for a model-based stabilizing gas lift con- 
3 0 troller for automatic and on-line control represents a 
number of inventive steps. First and foremost, our con- 
cept is able to stabilize the pressures, temperatures and 
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flow rates of a gas lifted well in an operating point 
that is unstable in open- loop (i.e., when no active con- 
trol is used) . The unstable production phenomena for a 
gas lifted well is eliminated, without increasing the 
5 mean gas lift injection rate, through active and continu- 
ous manipulation of the opening of the production choke 
and/or the opening of the gas injection choke as a dy- 
namic function of available process measurements. The 
model-based stabilizing gas lifted well controller pro- 
10 vides a way to stabilize gas lifted wells with different 
measurement devices (sensors) available for control pur- 
poses. 

The model based stabilizing gas lifted well controller is 
15 also characterized in that it allows that the control er- 
ror, which is a function of an externally given (optimal) 
reference operating point and the real operating point, 
at any time may be minimized with respect to a predefined 
model -based (integral) norm. 

20 

Brief Description of the Drawings 

Fig.l shows a schematic of a gas lifted oil well includ- 
ing the most common measurements . 

25 Fig. 2 shows unstable production versus time for a gas 
lifted oil well simulated in the multiphase simulator 
OLGA. 

Fig. 3 shows a typical lift performance relationship 
3 0 curves for both a stable and an unstable gas lifted oil 
well. 
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Fig. 4 shows a schematic of the invented control struc- 
ture . 

Fig. 5 shows an open loop simulation for an experimental 
5 gas lifted using a invented nonlinear gas lifted well 
model . 

Fig. 6 shows an illustration of input, output, and dis- 
turbances for a linear state-space gas lifted well model. 

10 

Fig. 7 shows a Bode plot of the transfer function from 
the gas injection choke to the wellhead pressure for an 
experimental gas lifted well. 

15 Fig. 8 shows a typical transfer function from the produc- 
tion curve to the casing pressure, both estimated value 
(x) and the fitted model (solid) . 

Fig. 9 shows a possible control structure for control of 
20 the wellhead pressure using the gas injection choke. 

Fig. 10 shows an illustration of input and output for the 
linear state-space model representing a locally linear 
stabilizing gas lifted well controller. 

25 

Fig. 11 shows the total invented model-based stabilizing 
gas lift controller consisting of several stabilizing gas 
lift controllers. 

3 0 Fig. 12 shows the structure for stabilization of the gas 
lifted well using measurements of wellhead pressure and a 
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model -based controller for manipulation of the gas injec- 
tion choke. 

Fig. 13 shows the structure for stabilization of the gas 
5 lifted well using measurements of bottom hole pressure 
and a model -based controller for manipulation of the gas 
injection choke. . 

Fig. 14 shows stabilization of WHP (left) and stabiliza- 
10 tion of the oil production rate (right) . 

Fig. 15 shows the structure for stabilization of the gas 
lifted well using measurements of LGR and a model -based 
controller for manipulation of the gas injection choke. 

15 

Fig. 16 shows stabilization of LGR (left) and stabiliza- 
tion of the oil production rate (right) . 

Fig. 17 shows the structure for stabilization of the gas 
2 0 lifted well using an estimate of LGR and a model -based 
controller for manipulation of the gas injection choke. 

Fig. 18 shows the structure for stabilization of the gas 
lifted well using an estimate of LGR and a model -based 

2 5 controller for manipulation of the gas injection choke. 

Fig. 19 shows the structure for stabilization of the gas 
lifted well using measurements of bottom hole pressure 
and a model -based controller for manipulation of the pro- 

3 0 duction choke. 
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Fig. 2 0 shows stabilization of bottom hole pressure 
(left) and stabilization of the oil production rate 
(right) . 

5 Fig. 21 shows the structure for stabilization of the gas 
lifted well using measurements of casing pressure and a 
model -based controller for manipulation of the production 
choke . 

10 Fig. 22 shows stabilization of casing pressure (left) and 
stabilization of the oil production rate (right) . 

Fig. 2 3 shows the multivariable control structure for 
stabilization of the gas lifted well using measurements 
15 of wellhead pressure, casing pressure, LGR and a model- 
based controller for manipulation of the production choke 
and the gas injection choke. 

Fig. 24 shows stabilization of casing pressure (left) and 
20 stabilization of LGR (right) . 

Fig. 25 shows stabilization of wellhead pressure (left) 
and stabilization of oil production rate (right) . 

2 5 Detailed Description of the Invention 

Referring now to Fig. 1, there is shown a schematic of a 
gas lifted oil well comprising a production tubing 1, a 
variable choke 2 for controlling the production rate and 
a variable choke 3 for controlling the injection of lift 

3 0 gas. The production choke 2 is mounted in a production 

tubing 1 and the gas injection choke is mounted on the 
gas injection flow-line 4. The well has either an annular 
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space 5 between the production tubing 1 and the casing 
10, or a separate gas supply tubing for injection of 
pressurized gas from the source 7 into the lower end of 
the production tubing l. One or more gas lift 
5 valves/nozzles 6 (unloading valves) are evenly distrib- 
uted alongside the production tubing 1; and during normal 
operation the gas is usually injected through the deepest 
of these valves. In the case of an annular space 5 (annu- 
lar space 5 is shown in Fig. 1) , a packer 8 is mounted 
10 near the lower end of the production tubing 1 to stop 

fluid from penetrating into the annular space 5 which is 
filled with pressurized gas. 

A wellhead assembly 9 makes the annular space leak- tight 
15 and may be equipped with devices for measurement of both 
wellhead pressure 11 and wellhead temperature 12 . In some 
cases it is also possible to measure the production flow 
rate 13 directly. In addition, the gas-conduit 4 may in- 
clude a device 14 for measurements of pressure in the an- 
2 0 nular space (casing/annulus pressure) , a temperature sen- 
sor for measurement of the annulus (casing) temperature 
18, and a device 15 for measurement of flow rate of the 
pressurized gas. Some wells are also equipped with sen- 
sors for measurements of the bottom hole pressure 16. At 
2 5 present time there exist no measurement devices for meas- 
urement of the gas injection rate 17 from annulus 5 to 
the production tubing 1 (referred to as Lift Gas Rate 
(LGR) ) , however, it is possible to estimate this quan- 
tity. 

30 

Fig. 2 shows a typical unstable production sequence (also 
referred to as heading, casing heading or heading limit 
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cycles) versus time from a gas lifted oil well. At the 
highest gas injection rates, the pressure drops in the 
tubing 1 is dominated by friction. If the gas-oil-ratio 
rises, the tubing pressure will increase which will re- 
5 duce the gas injection rate. This region therefore en- 
sures stable production. At low gas injection rates how- 
ever, the hydrostatic pressure gradient dominates the 
pressure drop in the production tubing 1. A small in- 
crease in gas-oil-ratio results then in a lower pressure 

10 in the tubing 1, which leads to a higher gas injection 

rate 17 from the annular space 5 into the production tub- 
ing 1 through the down-hole gas lift valve 6 . Since the 
gas rate is restricted by a gas injection choke 3 at the 
wellhead 9, the gas pressure in the annular space 14 will 

15 be reduced. After a time, the gas rate into the produc- 
tion tubing 17 will therefore be reduced, resulting in a 
lower oil production rate. At low gas injection rates the 
well is therefore intrinsically unstable in spite of the 
fact that wells are operated most efficiently on the up- 

2 0 ward slope of the lift performance relationship curve. 



Fig. 3 shows a typical lift performance relationship 
curves for gas lifted wells- The solid curve corresponds 
to injecting lift gas directly through the lowest gas in- 

25 jection valve 6 at constant rate. Obviously, this is not 
possible, and the dotted curve shows the resulting aver- 
age production when the lift gas is injected at constant 
rate through the gas injection choke 3. In both cases 
the resulting oil production rate are shown as a function 

3 0 of the amount of injected lift gas. Due to unstable pro- 
duction at low gas injection rates when injection the 
lift gas through the gas injection choke 3, it is seen 
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from Fig. 3 that the loss in average production is high 
for unstable production conditions (dotted curve) as com- 
pared to stable operating conditions (solid curve) . With 
the invented model-based stabilizing controller presented 
5 here it is claimed that the well in Fig. 3 may be oper- 
ated on or in infinitesimal distance from the stable 
solid curve. 

Fig. 4 shows a schematic of the invented model -based sta- 
10 bilizing gas lift controller structure. The model -based 
gas lift controller uses an externally given (optimal) 
reference point and one or more process measurements (or 
an estimate of these) to calculate the opening of the 
production choke 2 and/or the gas injection choke 3. The 
15 preferred mode for the externally given reference point 
is the specification of the optimal LGR 17 (gas rate 
through the down-hole gas injection valve) . 

In order to analyze and design controllers for gas lifted 
2 0 wells, our innovative idea is to design the stabilizing 
controller, and, if applicable, the estimator based on a 
dynamic model of the system. Therefore we have invented a 
structure for a simplified dynamic non-linear model based 
on physical principles of gas lifted wells suitable for 

2 5 controller and estimator design. The main purpose with 

this dynamic model is to describe the interactions be- 
tween the annular space 5 and tubing 1 which leads to the 
unstable behavior (heading limit cycles) at low and in- 
termediate gas injection rates. In addition it is neces- 

3 0 sary that the model becomes stable at high gas injection 

rates. 
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20 



25 



The idea is to use a simple model basically relying on 
three differential equations conserving mass in the tub- 
ing 1 and casing 5, and a couple of algebraic equations 
(of state) for approximating energy and impulse balances. 
At the cost of a more complicated, yet accurate, model, 
differential equations describing energy balances and im- 
pulse balances may also be included. 

The invented structure of a nonlinear dynamic gas lifted 
well model consists of: 

• Model of the pipes (casing 5 and tubing 1) : 

1. Three ordinary differential equations conserving 
masses in casing 5 and tubing -1. 

2 . Algebraic equations (of state) relating pressure, tem- 
perature, and liquid and gas holdup to each other in 
casing 5 and tubing 1. 

3 . Algebraic equations for pressure head. 

• Model of gas injection choke 3: An algebraic equation 
describing the relation between the pressure upstream 
and downstream the gas injection choke 3 and the mass 
flow rate through the choke. One possible equation to 
be used is : 



jection choke 3, w (3) is the opening of the gas injec- 
tion choke 3, and /? M(3) and /? J(3) are the upstream and 




Here, w i 



Ufigas 



is the mass flow rate through the gas in- 
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downstream pressures of the gas injection choke 3. C 0) 
is a constant parameter depending the gas injection 
choke used. 

• Model of the gas injection nozzle 6: An algebraic 
5 equation describing the relation between the pressure 

upstream and downstream the gas injection nozzle 6 and 
the mass flow rate through the nozzle 17. The equation 
will vary depending on the type of gas injection valve 
used. 

0 • Model of the production choke 2: An algebraic equation 
describing the relation between the pressure upstream 
and downstream the production choke 2 and the mass 
flow rate 13 of gas and liquid through the choke 2. 
One possible equation to be used is: 



W total = C (2) U (2) VP mix yjiPutf) PdJL2) ) 

Here, w total is the total mass flow rate through the pro- 
duction choke 2 , u (2) is the opening of the production 
0 choke 2, and /? w(2) and /? rf(2) are the upstream and down- 

stream pressures of the production choke 2 . C (2) is a 
constant parameter depending the production choke used 

The advantages with the invented model structure are : 

;5 

• It is compact (only a set of ordinary differential 
equations and algebraic equations) 

• It is able to capture the main dynamic behavior of gas 
lifted well both at low, medium (unstable operating 
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conditions) and high (stable operating conditions) gas 
injection rates. 

• The model may easily be linearized, meaning that it is 
suitable for controller and estimator design. 

5 • The parameters in the model can be tuned so that the 
model fits measured real time-series of pressures, 
temperatures, and flow rates from a gas lifted well, 

• The parameters in the model can be tuned so that the 
model fits simulated time-series of pressures, tem- 

10 peratures, and flow rates from a gas lifted well mod- 

eled in a rigorous multiphase simulator based on par- 
tial differential -algebraic equations 



Fig. 5 shows an open- loop simulation where the invented 
15 structure has been used to model an unstable experimental 
gas lifted well . The model has been implemented and simu- 
lated in MATLAB. As seen from the simulation results, the 
model is able to capture the oscillating limit cycle con- 
ditions, also known as casing heading, in the experimen- 
20 tal gas lifted well. 

A nonlinear dynamic gas lifted well model in accordance 
with the invented structure may be used directly as part 
of the model -based stabilizing gas lift controller shown 

25 in Fig. 4. However, it is sometimes difficult to design a 
model -based controller based on a nonlinear model. The 
preferred mode for utilizing the derived nonlinear model 
will therefore be linearization. To locally capture the 
dynamic behavior of an unstable operating point of a gas 

3 0 lifted oil well, the nonlinear model in accordance with 
the structure described above may be linearized in the 
current operating point of interest. Representing the lo- 
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cal dynamics of a gas lifted well using a linear state- * 
space model or, equivalently , a transfer function model 
then locally captures the dynamic behavior in the neigh- 
borhood of an unstable operating point. A linear state- 
5 space model of a gas lifted well will generally have the 
following format: 

jc(/) = Ax(t) + Bu(t) + Ed{t) 
y(t) = Cx(i) + Du(t) + Fd(t) 



10 Here, x(t) is the n x 1 state- space vector (in preferred mode 
n = 3 with the state-space representing mass of gas in the 
annular space 5 and tubing 1 and mass of liquid in the 
tubing 1) . u(t) is the 2x1 vector of controller inputs 
(opening of production choke 2 and gas injection choke 

15 3), d{t) is the kx\ vector of disturbances ( k = 5 in pre- 
ferred mode see description below) . Finally, y(t) is the 
pxl vector of outputs (in preferred mode ^corresponds 
to measured and/or estimated physical quantity) . The pa- 
rameters of the matrices B, C 9 D,E and F depend on the op- 

2 0 erating point where the nonlinear model has been line- 
arized. A more or less equivalent representation of this 
linear state -space model would be the corresponding 
transfer function representation: 



y{s) = G(s)u(s) 

G(s) = (C(sl - A)~ l B + D) 

Here, y(s) and u(s) are the Laplace transforms of y(t) and 
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The linear state-space gas lifted well models have the 
following inputs: 

• Opening of the gas injection choke 3 
5 and /or 

• Opening of the production choke 2 

and it has the following outputs: 

10 • Wellhead pressure 11 

and/or 

• Bottom hole pressure 16 
and/or 

• Casing pressure 14 
15 and/or 

• Mass rate of gas through gas injection valve 17 
and /or 

• Casing temperature 18 
and/ or 

2 0 • Mass rate of gas through gas injection choke 15 

and the following disturbances: 

• Pressure and temperature upstream the gas injec- 
25 tion choke 3 

and/or 

• Pressure and temperature in the reservoir 
and/or 

• Pressure downstream the production choke 2 

30 
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Fig. 6 shows all possible inputs, outputs and distur- 
bances for the linear state-space gas lifted well model. 

We have invented two ways of generating these kind of 
5 linear gas lifted well models: 

1. As already alluded to, by numerical or analytical 
linearization of a nonlinear dynamic gas lifted well 
model in accordance with our invented structure de- 
10 scribed above. 

Example ; From numerical linearization of the nonlinear 
gas lifted well model simulated in Fig. 5, we find 
that a typical transfer function from e.g. the gas in- 
15 jection choke 3 to the wellhead pressure 11 in an un- 

stable operation point is. given by: 



58.2 (& + 6.3 10~ 4 ) 

G(s) - 



(s-8.8-10- 4 )(s-2.9.10 3 )(s-h73.10- 2 ) * 

20 It is immediately seen that there are two poles in the 

right half plane, which means that the system is un- 
stable in open loop. Bode plots of G(s) are shown in 
Fig . 7 . 

25 2 . By closed-loop identification experiments on a gas 

lifted oil well modeled in a multiphase pipeline simu- 
lator (OLGA) where the closed- loop system is stable in 
the operating point in question. 

3 0 Example: By using ABB'S MATLAB / OLGA link we have in- 

vented a way to run such experiments from MATLAB. An 
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5 



example of a transfer function identified in this way 
is given in Fig. 8, where the transfer function from 
the production choke 2 to casing head pressure in the 
nominal operation point of a typical gas lifted oil 
well have been identified using this approach. 



Used in combination with advanced techniques from control 
theory, the linear local gas lifted well models (as de- 
scribed above) can be used to design model -based linear 
10 locally stabilizing gas lift controllers. In this way, an 
(optimal) operating point that is unstable in open- loop 
(i.e. without active control), becomes locally stable in 
closed-loop (i.e. when the stabilizing gas lift control- 
ler is actively used) . 

15 

Fig. 9 shows an example where a simple SISO (Single Input 
Single Output) control structure is used for controlling 
the wellhead pressure 11 using the gas injection choke 3. 
Based on the typical transfer function from the gas in- 
2 0 jection choke to the wellhead pressure 11, it is 

straightforward to design a stabilizing controller. In 
this example, a controller given e.g. by: 



25 Resulting in the following stable closed-loop poles: 



Controllers) = 5.70e - 5(1 + 



1 



) 



3600s 



0.034 + 0.042i 



0.034 - 



0 . 042i 



0. 00083 



30 



0.00023 
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Generally, our linear locally stabilizing gas lift well 
controllers will become more complex than the simple 
SISO-controller illustrated above, and therefore we rep- 
resent them in linear state- space form, or, equivalently , 
5 as transfer functions. A linear state-space model of a 

linear stabilizing gas lifted well controller will gener- 
ally have the following format: 

x c (0 = A c x(t) + B c y{t) + E c r(t) 
u(i) = C c x c (/) + D c y(t) + F c r(t) 

10 

Here, x c (t) is the wxl controller state -space vector, u(t) 
is the 2x1 vector of controller outputs (opening of pro- 
duction choke 2 and gas injection choke 3) . y{t) is the 
px\ vector of controller inputs (in preferred mode y(t) 

15 corresponds to measured and/or estimated physical quan- 
tity) . /*(?) is the qx\ vector of the externally given op- 
erating point. The parameters of the matrices 
A c ,B c ,C c ,D C9 E c a.n<l F c depend on the parameters in the lin- 
ear state-space gas lifted well model in the operating 

20 point in question. These linear state-space models, rep- 
resenting the linear model-based linear stabilizing gas 
lifted well controllers, have the following inputs: 

• Wellhead pressure 11 

2 5 and/ or 

• Bottom hole pressure 16 
and /or 

• Casing pressure 14 
and/or 

3 0 • Mass rate of gas through gas injection valve 17 
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and/or 

• Casing temperature 18 
and /or 

• Mass rate of gas through gas injection choke 15 
5 and/or 

• Pressure and temperature upstream the gas injec- 
tion choke 

and/or 

• Pressure and temperature in the reservoir 
10 and/or 

• Pressure downstream the production choke 
and 

• Externally given (optimal) reference operating 
point 



15 



and the following outputs: 



• Opening of the gas injection choke 3 
and /or 

2 0 • Opening of the production choke 2 

Fig. 10 shows all possible inputs and outputs for the 
linear state- space models representing a linear stabiliz- 
ing gas lifted well controller. 

25 

In order to generate globally model-based stabilizing gas 
lifted well controllers, our innovative idea is to com- 
bine the model -based linear locally stabilizing gas 
lifted well controllers described above. Each total 
30 model-based gas lifted well controller will then consist 
of a family of model -based linear stabilizing control - 
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lers, each of which will be valid in a predefined neigh- 
borhood of an open- loop unstable operating point, and 
switching between the controllers will occur based on 
predefined logical rules. In addition to logical switch- 
5 ing rules, we also include logic to prevent integrator 
windup whenever integral action is included in the con- 
troller. Fig. 11 illustrates the total gas lifted well 
controller, including several linear model-based stabi- 
lizing controllers . 

10 

Examples 

Several control structures, in line with our general con- 
cept shown Fig. 4, with the ability to stabilize a gas 
lifted well in an operating point that is unstable in 
15 open- loop, are described in the examples that now fol- 
lows. What is common for all these model-based controller 
concepts is that the heading phenomena is eliminated 
through active and continuous manipulation of the opening 
of the production choke 2 and/or the opening of the gas 

2 0 injection choke 3. 

In the examples, the gas lifted well is modeled in the 
multiphase simulator OLGA and the model-based gas lift 
controller is implemented in MATLAB „ The experiments have 
25 been done in MATLAB using ABBs MATLAB/OLGA link. 

Example 1 

Using only measurements of pressure in the production 
tubing 1 as input to a model based gas lift controller, 

3 0 the gas lifted well may be stabilized only through dy- 

namic manipulation of the gas injection choke 3. Pressure 
in production tubing 1 may be measured anywhere between 
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the bottom of the well, i.e. bottom hole pressure 16, to 
the wellhead 9, i.e. wellhead pressure 11. Fig. 12 and 
Fig. 13 shows the controller structure using measurements 
of wellhead pressure 11 and bottom hole pressure 16. 

5 

Simulations are performed where the controller manipu- 
lates the gas injection choke 3 in order to stabilize the 
wellhead pressure 11. Fig. 14 shows simulation results 
for this control concept. The controller starts after 
10 eight hours of open-loop simulations. 

Example 2 

Using only measurements of the LGR 17 as input to a model 
based gas lift controller, the gas lifted well may be 
15 stabilized only through dynamic manipulation of the gas 
injection choke 3. The controller structure using this 
measurement is shown in Fig. 15. 

Results from simulations where the model -based controller 

2 0 is used for manipulation of the gas injection choke 3 for 

stabilization of LGR 17, are shown in Fig. 16. The con- 
troller starts after eight hours of open- loop simula- 
tions . 

25 Example 3 

If measurements of LGR 17 through the active gas injec- 
tion valve 6 are not available, a non- linear model -based 
dynamic estimator may be used to estimate this rate. The 
estimator may use the measurements of the gas injection 

3 0 rate through the gas injection choke 15, temperature in 

casing 18 and pressure in casing 14. In addition to these 
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measurements, the estimator may use the opening of the 
gas injection choke 3 itself. 

Using only an estimate of LGR 17 (based upon the non- 
linear model -based dynamic estimator described above) as 
input to a model -based dynamic gas lift controller, the 
gas lifted well will be stabilized only through dynamic 
manipulation of the gas injection choke 3. The lift gas 
rate 17 is controlled indirectly when using the estima- 
tor. The controller structure using the estimator is 
shown in Fig. 17. 

Using measurements of the pressure in the production tub- 
ing 16, 11 the pressure in the casing 14 and the opening 
of the gas injection choke, LGR 17 may be estimated. 
Based upon this estimate, LGR 17 may be controlled indi- 
rectly only through dynamic manipulation of the gas in- 
jection choke 3. The controller structure using an esti- 
mate of lift gas rate 17 is shown in Fig. 18. 

Example 4 

Using only measurements of bottom hole pressure 16 as in- 
put to a model based gas lift controller, the gas lifted 
well will be stabilized only through dynamic manipulation 
of the production choke 2. The controller structure using 
this measurement is shown in Fig. 19. 

Results from simulations where the model -based controller 
is used for manipulation of the production choke 2 for 
stabilization of the bottom hole pressure 16 is shown in 
Fig. 20. The controller starts after three hours of open- 
loop simulations. 
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Example 5 

Using only measurements of pressure in casing 14 as input 
to a model based gas lift controller, the gas lifted well 
5 will be stabilized only through dynamic manipulation of 
the production choke 2. The controller structure using 
this measurement is shown in Fig. 21. 

Results from simulations where the model -based controller 
10 is used for manipulation of the production choke 2 for 

stabilization of the bottom hole pressure 16 is shown in 
Fig. 22. The controller starts after three hours of open- 
loop simulations. 

15 Example 6 

Several measurements, and/or an estimate of one or more 
of these, may be used as input to a multivariable model - 
based gas lift controller. On possible structure for this 
multivariable controller is shown in the example below. 

20 

Using measurements of pressure in the casing 14, pressure 
at the wellhead 11 and LGR 17 as input to a multivariable 
model-based gas lift controller, the gas lifted well will 
be stabilized through dynamic manipulation of both the 
25 production choke 2 and the gas injection choke 3. The 

controller structure using these measurements is shown in 
Fig. 23. 

Results from simulations where the model -based controller 
3 0 is used for manipulation of the production choke 2 and 
the gas injection choke 3 for stabilization of the gas 
lifted well, is shown in Fig. 24 and in Fig. 25. The 
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multivariable controller starts after fourteen hours and 
after sixteen hours the setpoint for LGR 17. is ramped 
from 0.6 kg/s to 0.8 kg/s. 

5 Example 7 

From closed- loop experiments, we have invented a way to 
tune the controller parameters on-line (c.f. page 15). 
Tests have successfully been performed for on-line tun- 
ing . 

10 . 

To determine the optimum setpoint value, a logic sequence 
combined with a stepwise approach has been used. Using 
this method to determine the optimum setpoint value, no 
other inputs will be required from the well. 

15 
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Claims 

1. Method for controlling the production flow rate of 
an oil well, said well comprising a production tubing (1) 
5 with at least one production choke (2) and gas injection 
means including at least one gas injection choke (3) , one 
or more of the chokes being controlled as a function of 
process measurements , 

characterized in that 
10 - the well, the production of the well and the conditions 
in the gas injection means are monitored, measured and/or 
estimated, 

- pressures, temperatures and flow rates are actively 
stabilized even in an open- loop unstable operating point 

15 by means of a model -based control system comprising a 
stabilizing controller, and 

- said choke or chokes is continuously manipulated by the 
stabilizing controller through feedback control as a dy- 
namic function of model-based estimates of process pres- 

2 0 sures, temperatures and/or flow rates, and/or said proc- 

ess measurements. 

2 . Method according to claim 1 , 

characterized in that a mathematical dy- 
25 namic model is made of the well system, the model being 
comprised in the model-based control system in associa- 
tion with the stabilizing controller and having the abil- 
ity to describe and recreate the unstable limit cycles 
that may occur in pressures, temperatures and flow rates 

3 0 in the production tubing (1) and/or gas supply means (5) 

included in the gas injection means for supply of pres- 
surized gas to the lower end of the production tubing. 
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3 . Method according to claim 2 , 

characterized in that the stabilizing 
controller is designed and tuned based on the model. 

5 

4. Method according to claim 2, 

characterized in that the mathematical 
dynamic model of the well system is non-linear, in order 
to capture the behavior over a wide operating range, and 
10 based on ordinary differential and algebraic equations. 

5. Method according to claim 2 or 4 , 

characterized in that one or more of the 
parameters in the model are adjusted in order to fit the 
15 model to measured time series of pressure, temperature 
and flow rates from a well. 

6. Method according to claim 2 or 4, 

characterized in that one or more of the 
20 parameters in the model are adjusted in order to fit the 
model to simulated time series of pressures, temperatures 
and flow rates from a well which is modeled in a rigorous 
multiphase pipeline simulator based on partial differen- 
tial-algebraic equations. 

25 

7 . Method according to claim 2 , 

characterized in that the model is a com- 
bination of a number of linear state -space 
models, each linear state-space model being represented 
3 0 by a set of system matrices or an equivalent representa- 
tion, each linear state-space model simulating the dy- 
namic behavior of an oil well in the neighborhood of an 

BNSDOCID: <WO 000071 5A1_I_> 



WO 00/00715 PCT/NO99/00217 

26 

open- loop unstable operational point, each linear state- 
space model comprising one or both of the following in- 
puts : 

-opening of the gas injection choke (3) , 
5 -opening of the production choke (2) , 

and comprising one or more of the following outputs: 
-wellhead pressure (11) , 
-bottom hole pressure (16) , 

-casing pressure (14) /pressure in gas supply tubing, 
10 -mass rate of gas through gas injection valve (17) , 

-casing temperature ( 18 ) /temperature in gas supply tub- 
ing, 

-mass rate of gas through gas injection choke (15) , 
and, if necessary, one or more of the following distur- 
15 bances : 

-pressure and temperature upstream the gas injection 
choke , 

-pressure and temperature in the reservoir, 
-pressure downstream the production choke. 

20 

8 . Method according to claim 7 , 

characterized in that each linear model 
is derived through a numeric or algebraic linearization 
of a non- linear dynamic model of the well system, with 
25 the ability to capture the behavior over a wide operating 
range and being based on ordinary differential and alge- 
braic equations. 

9. Method according to claim 2 and 7, 

30 characterized in that each linear state- 
space model is identified through experimental closed- 
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loop perturbation of a well system which is modeled in a 
multiphase pipeline simulator. 

10- Method according to claim 2, 
5 characterized in that the stabilizing 

controller is represented as a combination of a number of 
linear state-space models, each linear state-space model 
being represented by a set of system matrices or an 
equivalent representation, each linear state- space model 

10 simulating the dynamic behavior of a linear stabilizing 
well controller in such a way that an open- loop unstable 
operating point for pressures, temperatures and flow 
rates is stabilized in closed-loop in the neighborhood in 
which the linear state-space model is valid, each linear 

15 state- space model comprising one or more of the following 
inputs: 

-wellhead pressure (11) , 
-bottom hole pressure (16) , 

-casing pressure (14) /pressure in gas supply tubing, 
20 -mass rate of gas through gas injection valve (17), 

-casing temperature (18) /temperature in gas supply tub- 
ing, 

-mass rate of gas through gas injection choke (15) , 
and comprising one or more of the following outputs: 
25 -opening of the gas injection choke (3), 
-opening of the production choke (2) . 

11. Method according to claim 7 and 10, 

characterized in that the linear state- 
3 0 space models comprising the stabilizing controllers are 

derived based on the linear state-space models comprising 
the dynamic well model . 
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12. Method according to claim 8 or 9, 

characterized in that the linear state- 
space models comprising the stabilizing controllers are 
5 derived based on the linear state -space models comprising 
the dynamic well model. 

13. Method according to claim 2, 

characterized in that the stabilizing 
10 comtroller is represented by a set of non-linear ordinary 
differential equations and/or algebraic equations in or- 
der to stabilize the well system over a wide operating 
range . 

15 14. Stabilizing well controller for controlling the pro- 
duction flow rate of an oil well, said well comprising a 
production tubing (1) with at least one production choke 
(2) and gas injection means including at least one gas 
injection choke (3) , the chokes being controlled as a 

2 0 function of process measurements and/or model -based esti- 
mates of pressures, temperatures and flow- rates, 
characterized in that said stabilizing 
controller is comprised within a model -based control sys- 
tem which further comprises a dynamic mathematical model 

2 5 of the well system and that said stabilizing controller 

and mathematical model are associated. 

15. Stabilizing well controller according to claim 14, 
characterized in that said stabilizing 

3 0 controller is arranged to be designed and tuned based on 

the mathematical model and/or process measurements. 
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16. Stabilizing well controller according to claim 14 or 
15, 

characterized in that said stabilizing 
controller comprises a number of stabilizing controllers, 
5 each of which being valid in a predefined neighborhood of 
an open- loop unstable operating point, and that the con- 
troller comprises or is associated with means for switch- 
ing between said controllers based on predefined logical 
rules comprised in the mathematical model . 

10 

17. Stabilizing well controller according to claim 14, 
15 or 16, 

characterized in that it comprises built- 
in logic and/or non-linearities to prevent integrator 
15 windup and input saturation. 

18. Stabilizing well controller according to claim 14 or 
17, 

characterized in that it manipulates the 
20 opening of the gas injection choke (3) using the measure- 
ment of the pressure in the production tubing (1) as in- 
put . 

19. Stabilizing well controller according to claim 18, 
25 characterized in that it uses the meas- 
urement of the bottom hole pressure (16) as input. 

20. Stabilizing well controller according to claim 18, 
characterized in that it uses the meas- 

30 urement of the wellhead pressure (11) as input. 
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21. Stabilizing well controller according to claim 14 or 
17, 

characterized in that it manipulates the 
opening of the gas injection choke (3) using the measure - 
5 ment of the lift gas rate from casing (17) /gas supply 
tubing to the production tubing as input . 

22 . Stabilizing gas lift controller according to claim 
14 or 17, 

10 characterized in by a non- linear dynamic 
well measurement filter (model-based estimator) , said es- 
timator being arranged to utilize the controlled measure- 
ments of the gas injection rate through the gas injection 
choke (15) , temperature and pressure in the casing (18, 

15 14) /gas supply tubing, and said estimator estimating the 
rate of lift gas through the active gas injection valve 
(17). 

23 . Stabilizing well controller according to claim 14 or 
20 17 , 

characterized in that said controller 
based on an estimate from a non- linear gas lift filter is 
arranged to manipulate the opening of the gas injection 
choke (3) in order to indirectly control the lift gas 
25 rate from casing/gas supply tubing to the production tub- 
ing (17) . 

24. Stabilizing well controller according to claim 14 
and 17, 

30 characterized in that said controller is 
arranged to, based on an estimate based on measurements 
of pressure in the production tubing (1) and pressure in 
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casing (14) /gas supply tubing, manipulate the opening of 
the gas injection choke in order to indirectly control 
the lift gas rate from casing/gas supply tubing to the 
production tubing (17) . 

5 

25. Stabilizing well controller according to claim 14 or 
17, 

characterized in that said controller is 
arranged to, based on the measurements of the bottom hole 
10 pressure (16) as input, manipulate the opening of the 
production choke (2) . 

26. Stabilizing well controller according to claim 14 or 
17, 

15 characterized in that said controller is 
arranged to, based on the measurements of pressure in 
casing (14) /gas supply tubing as input, manipulate- the 
opening of the production choke (2) . 

20 27. Stabilizing well controller according to claim 14 or 
17 > . 

characterized in that said controller is 
arranged to, based on the measurements of pressure in 
casing (14) /gas supply tubing and at the wellhead (11) as 
25 input, manipulate the opening of both the production 
choke (2) and the gas injection choke (3) . 

28. Stabilizing well controller according to any of the 
previous claims 14 - 27, 
30 characterized in that said controller is 
arranged to at any time minimize the deviation between 
the optimal reference operating point and the real oper- 
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ating point (control error) , with respect to a given time 
horizon . 



29. Stabilizing well controller according to claim 28, 
5 characterized in that said controller is 
arranged to give the optimal reference operating point at 
an optimal gas injection rate from casing/gas supply tub- 
ing to the production tubing (17) . 

10 30. Stabilizing well controller according to any of the 
previous claims 14-27, 

characterized in that said controller is 
arranged to adjust the parameters in the controller on- 
line through closed- loop perturbations. 

15 



BNSDOCID: <WO 0O0O7 1 5A 1 _l_> 



WO 00/00715 



1 / 25 



PCT/NO99/00217 



12 





.13 





i 



QAA-' 

14 18 



.17 



BNSDOCID: <WO 00007 15A1J_> 



WO 00/00715 



PCT/NO99/00217 




o 
o 



OX) 



CO 



CM 



O 
CM 



CO 



CNJ 



00 



(s/B>j) e;ej uoi^onpojd ijq 



BNSDOCID: <WO 000071 5A1_L> 



WO 00/00715 



.3 / 25 



PCT/NO99/00217 




WO 00/00715 



PCT/NO99/00217 



4 / 25 



CO 

co 



i2 

CD 

E 

0> CD 
O b- 

t: co 
CL co 

CD 



CO = 
CD O 

o 



E g 

la 




BNSDOCID: <WO 00007 15A1_L> 



WO 00/00715 



5 / 25 



PCT/NO99/00217 




WO 00/00715 



PCT/NO99/00217 



6/25 



CL 
I— 



r 



T3 
Q> 



CD 



CL 
_CD 
O 
XT 

E 
o 

o 

CD 



CL 

.£ 

CO 

co 
O 



IT 
O 



Cl 

E 
© 

cr 

CO 

ro 
O 



© 

o 
x: 
O 

c 
o 

o 



CO 
CO 

O 

x: 
cr 

2 



CD 

a: 



CO 



CD 
O 

« 8. • 

o co 3* 
c » o 
.r o> ^ 

05 



© 

o 

-C 

O 
c 
.2 
o 
3 

o 



o 
cr 
jo 
c 

Q> 

o 



© 

O 
-C 

O 

c 
.o 

o 

CD 



CO 
CD 
O 

o 

o> 
c 

© 

CL 

O 



CD 
O 

x: 
O 



CO 

O 
E 

CO 



E 
c 

CD 
CO 

£ 

CL 



i 

CD 
CO 
CD 

a: 



o. 
E 

CD 
h~ 
TD 
C 
CD 



3 
CO 
CO 

s? 

CL 



CD 

O 

XT 

O 
c 
o 



73 

e 

CL 

E 

CO 

CD 

To 
c 

o 
Q 

3 

CO 
CO 

CD 



3 
CL 



fa 



C/) 

o 

< 
CD 

cr 

cr 
i — 

CO 

o 



BNSDOCID: <WO 000071 5A1_L> 



WO 00/00715 



7/25 



PCT/NO99/00217 




WO 00/007 1 5 PCT/NO99/002 1 7 

8 / 25 




BNSDOCID: <WO 000071 5A1 I > 



WO 00/00715 



PCT/NO99/00217 



9 / 25 



co 

CO 

a? 

CL 

T5 
CO 
CD 



CD 



• MM 



CD 

O 

o -c 
0)0 

E c 
c .2 

°f 

CO 
CT3 

o 




. CO 
o CO 

C CL 

O "D 

Q. CO 

O © 
00 = 
CD 



BNSDOCID: <WO 000071 5A1_I_> 



WO 00/00715 



PCI7NO99/00217 



10 / 25 



Q. 

o 




O 

J 13 S 

CO 



2 

3 

CO 

co 



CD 

0) 

XT 

a3 



CD 



i— 
3 

CO 
CO 

Q_ 

JD 

O 
x: 

E 
o 

IS 

o 

CD 



C 
CO 
CD 

O 



or 
o 



co 



a 
E 

a 
c 
"to 
CO 



^~ 
O 

o 
x: 
O 



o 

JE 

CO 
CD 

o 

x: 
en 

3 
O 

CD 

q: 

CO 
CO 
CD 



CUD 



BNSDOCID: <WO 000071 5A1_L> 



WO 00/00715 PCT/NO99/00217 

11/25 




WO 00/00715 



PCT/NO99/00217 



12 / 25 



O <p 
C 3 

to t: 
co co 

a? 




CO 

w- O 
O JZ 
D)0 

c c 
c o 
a> = 
CL <-> 
O ^ 

CO 

O 




0X3 



Q> 

w 

to o 

w c 

a> cd 
*- 

O. 

-o "a? 

co a: 

a> 
5 



BNSDOC1D: <WO 000071 5A1_L> 



WO 00/00715 



PCT/NO99/00217 



13 / 25 



fa 



o 2? 

=> 

C M 

£ 2! 

0) CL 

« 1 

o 
o 

CD 




CO 
CD 

= c 
c o 

CD S3 

a. o 
O ^ 



w 

CD 

a 




a> 



<o CP 

to o 

a> c 

. >* 0) 



CD 



o 

CD 



BNSDOCID: <WO 000071 5A1_1_> 



WO 00/00715 



PC77NO99/00217 




WO 00/00715 



PCT/NO99/00217 



15 / 25 



So: 

ro — i 

0) 

2 




CO 

w- O 

O -C 

o>U 
£ c 
c o 

2-1 
of 

CO 

ca 
O 




OJD 



o 

—J 

1— 

o 
c 

0T 



BNSDOCID: <WO 000071 5A1J_> 



WO 00/00715 



16 / 25 



PCT/NO99/00217 




BNSDOCID: <WO 00007 1 5A 1 _ L> 



WO 00/00715 



PCT/NO99/00217 



17 / 25 




• * 



WD 

fa 



O 



0) 

o 
c 



BNSDOCID: <WO 000071 5A1_I_> 



WO 00/00715 



PCI7NO99/00217 




WD 




cr 
o 
_j 

05 

O 

c 

cr 



BNSDOCID: <WO 000071 5A1J_> 



WO 00/00715* 



PCT/NO99/00217 



19 / 25 



a> to 
E £ 
2 cl 

co 2 

fj 

o 
CD 




o o 

jf 6 

Q> O 

ol 
e 



si 




3 

CO Q> 

tr> o 
© C 

a £ 
© 3 

IS 

£ • 
o 

o 

CO 



BNSDOCID: <WO 00007 15A1_I_> 



WO 00/00715 PCT/NO99/00217 

20 / 25 



Fig. 20: 
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